Thermally activated conductivity in gapped bilayer graphene 
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FIG. 1: Panel (a) shows bilayer graphene's crystal structure 
and which-layer pseudospin orientation. Panel (b) shows the 
lowest two bands and band gap size A. The chemical potential 
fj. is counted from the middle of the band gap. 
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This is a theoretical study of electron transport in gated bilayer graphene — a novel semicon- 
ducting material with a tunable band gap. It is shown that the which-layer pseudospin coherence 
enhances the subgap conductivity and facilitates the thermally activated transport. The mechanism 
proposed can also lead to the non-monotonic conductivity vs. temperature dependence at a band 
gap size of the order of lOmeV. The effect can be observed in gapped bilayer graphene sandwiched 
in boron nitride where the electron-hole puddles and fiexural phonons are strongly suppressed. 



I. INTRODUCTION 

Graphene^ is often considered' as a most promis- 
ing material for future semiconductor industry. In- 
deed, it demonstrates high carrier mobility even at room 
temperature^ and is suitable for mass production thanks 
to the chemical vapor deposition technique developed 
recentlyi^i^ However, pristine graphene^ does not have 
a band gap which is a crucial ingredient for the field ef- 
fect transistor functionality. It is possible to open the 
gap in bilayer graphene by applying an external electric 
field perpendicular to the sample, see fig. [T] The effect 
was predicted by McCann^ and experimentally proven in 
ref;^. Note that it is also possible to open a gap between 
hole and conduction bands in bilayer graphene by means 
of an appropriate chemical dopingi^ 

In order to control the band gap and carrier density 
independently the double-gated graphene devices have 
been utilized^rJ^. The most striking feature observed is 
that the band gap obtained by infrared spectroscopyi^i^^ 
turns out to be much too large to fit the thermally ac- 
tivated conductivity measurements. There are a few at- 
tempts to resolve this discrepancy. An earlier modeli^ 
suggests the formation of midgap states in which charge 
carriers are localized. The band edge moves locally fur- 
ther into the gap and a hopping mechanism dominates 
the conductioni^iiS The most recent approachi^ employs 
fiuctuations of the charged impurity potential separat- 
ing the electron and hole puddles. Indeed, the first ex- 
perimental observations^"'^ ^ of the insulating behavior in 
gapped bilayer graphene have been made in the devices 
with graphene flakes placed directly on the Si02 sub- 
strate. The substrate impurities are known to cause siz- 
able potential fluctuations which lead to the formation 
of electron-hole puddles at low carrier densities JJ" If the 
substrate potential fluctuations are strong enough then 
the small effective band gap is expected to be due to 
the percolation through the charge inhomogeneities over- 
whelming the real spectral gap. The relevance of this 
mechanism to the subgap conductivity is unquestionable 
as long as graphene is placed on the Si02 substrateJ^ In 
recent experimentsiS- carried out on suspended double- 
gated bilayer graphene the electron-hole puddles are ex- 
pected to be suppressed; nevertheless, the activation en- 
ergy deduced from the transport measurements is still 



smaller than the band gap size. An alternative modeli^ 
suggests that the edge transport plays an important role 
in these measurements.^^ The phenomenon originates 
from non-trivial topological properties of the electronic 
band structure in graphene which are similar to those in 
spin-orbit induced topological insulators;^ 



A question addressed in this paper is whether there is 
another mechanism responsible for the substantial sub- 
gap conductivity which can manifest itself in gapped bi- 
layer graphene sandwiched in boron nitride i^^i^^ Such 
graphene samples are practically insusceptible to the 
environment making the substrate much less impor- 
tant. Moreover, the electron-hole puddles can be com- 
pletely screened out in double-layer systems similar to 
those recently reported in refi^. The edge transport, 
if any, can be precluded in Corbino geometry which 
has been already utilized in recent experiments carried 
out on double-gated bilayer graphenei^ Using the pseu- 
dospin coherence concept we predict that the subgap con- 
ductivity contribution does not vanish completely even 
though all abovementioned mechanisms are excluded, see 
flgs. I3I4I The signature of the mechanism in question 
is the non-monotonic conductivity vs. temperature de- 
pendence at a band gap size of a few tens of meV, see 
flgs. I5I6I This non-monotonic dependence could not be 
explained within conventional model-i^ii^ where disor- 
der renormalizes the band gap to a smaller value just by 
locally raising or lowering the band edges. 
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These are not the valence 
band eigenstates. 




FIG. 2: Panel (a) schematically shows the in-plane pseudospin 
component oc sini9fe for a valence band electron while it is 
moving through disordered gapped bilayer graphene. Each 
scattering event can be seen as a classical "measurement" 
which changes the carrier's wave function to a valence band 
eigenstate |ak) in accordance with the momentum dependent 
pseudospin texture shown on the right. At the same time the 
direction of particle's motion is changed and since the older 
pseudospin orientation a does not correspond to its eigenvalue 
with the new momentum k' the particle gets out of its eigen- 
state immediately after the scattering event. The subsequent 
collision returns the particle to its eigenstates with k' and a' 
but alters the direction of motion again. Thus, the valence 
electron is obviously not in its equilibrium valence eigenstate 
while it is moving between the scatterers and, therefore, can 
make a contribution to the total conduction even though the 
valence band is fully occupied. Panel (b) corresponds to the 
valence band electrons described by eq. ([2} where the pseu- 
dospin orientation does not depend on the direction of motion. 
The particle remains in its eigenstate after each scattering 
event. The valence band eigenstates can not conduct here as 
long as the valence band is fully occupied. 



II. CONCEPT 

We argue that the difference between effective (trans- 
port) and actual (spectral) gaps is an intrinsic property 
of gapped bilayer graphene following from the minimal 
two-band effective Hamiltonian already employed in ref;^. 
The Hamiltonian can be written as ^ hk ■ (J, where 
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(i cos 2ifi + y sin 2(p) + zU, 



(1) 



and a are Pauli matrices representing the pseudospir^ 
degree of freedom for carriers in bilayer graphene 
which originates from its peculiar crystal lattice shown 
in fig. [IJa) with the cr^ -pseudospin projection refer- 
ring to the layer index. Here, m is the effective 



mass, k is the two-component particle momentum, 
tan(/3 — ky/kj., and A = 2U is the band gap. The 

eigenvalues of Hq are iJ^k — n\J (^^^)^ + U'^ with 
K = ± being the band index, and the eigenstates 
are VKk(r) = e'^^^lx^k) with the spinors \x+k) = 
(cos^,sin^e2*^)^, \x-k) = (sin ^, - cos ^e^'^)^, 

where cos^k = U / ^JJ^^^f~+lP . The bands E^u are 
shown in fig.HJb). In order to pinpoint the mechanism re- 
sponsible for the transport gap renormalization we com- 
pare the pseudospin-momentum coupled model ([T]) with 
the decoupled one in which Hq = h'^. ■ a, where 
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zU. 
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The two models have the same energy spectrum but the 
the eigenstate spinors do not depend on the direction of 
particle's motion here and read |x+fc) = (cos sin ^)"^, 
\x'-k) — (sin — cos ^)-^. This only difference between 
two model Hamiltonians (H]) and Q leads to the drastic 
change in the subgap conductivity behavior: The subgap 
conductivity of graphene does not vanish even at zero 
temperature, see fig. [3l whereas it docs so within decou- 
pled band model, as shown in fig. 31 

The mechanism can be understood from fig. [21 Due 
to the pseudospin-momentum coupling in graphene the 
particle necessarily gets out of its valence band eigenstate 
while moving between two subsequent collisions with the 
scatterers. The resulting wave function, to a certain ex- 
tent, can be seen as a superposition between valence and 
conduction band states. (The conduction band state ob- 
viously represents an evanescent wave function as long as 
the energy is below the bottom of conduction bandit) 
As consequence the electron and hole states become en- 
tangled and in that way can facilitate the conductivity 
making the effective band gap smaller than the actual 
one, see figs. 

Note that such an entanglement has nothing to do with 
the electron- hole pairs. The electron- hole pairs are en- 
tirely classical objects and occur in both graphene and 
conventional semiconductor material as soon as the tem- 
perature reaches the level high enough to excite the va- 
lence electrons across the band gap. The interband en- 
tanglement is certainly of quantum mechanical nature. 
This phenomenon, as many other effects related to the 
quantum mechanical coherence, is sensitive to temper- 
ature. In some cases one can observe the competition 
between the temperature-dependent pseudospin decoher- 
ence and thermal activation of the electron-hole pairs 
which results in the non-monotonic conductivity vs. tem- 
perature dependence, see figs. jSHHl 



III. METHODS 

To evaluate the dc conductivity a we follow the proce- 
dure described i n^^i^^ and start from the finite-size Kubo 
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FIG. 3: Zero-temperature conductivity (for given spin and 
valley) of gapped bilayer graphene (band gap size A = 
87meV) in presence of the 5-correlated disorder with the 
strength uq = 2.74 • 10~^* eVcm'^. The concentration of scat- 
terers = Ns/L^ (with L = 1.8 x 10~^ cm being the sam- 
ple size) is different for each curve: (a) 0.54 ■ lO^^cm-2, (b) 
0.81-10^^ cm"^, (c) 1.62-10^^ cm~^. The coupling t? is chosen 
to be equal to 1Q6E, where 5E = 2i\f? jlT'm. The inset shows 
that the dependence of both metallic and subband conductiv- 
ities on r\ is relatively weak in this case. 
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FIG. 4: Zero-temperature conductivities for gapped bilayer 
graphene (solid lines) and conventional intrinsic semiconduc- 
tor (dotted lines) described by eqs. ([TJ and ([2} respectively. 
Note that the latter drops down to zero as soon as the Fermi 
energy level reaches the bottom of the conduction band. Dis- 
order parameters are the same as in fig. [3] for curve (b). The 
inset shows how the conductivity curve changes under scal- 
ing. The disorder concentration is chosen to be the same 
for all L's considered. One can see that the conductivity gets 
less sensitive to scaling for larger L's considered in this work. 
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where is the finite-size system area, 77 is the coupling 
to source and drain reservoirs, v is the velocity oper- 
ator, is the Fermi-Dirac distribution function, and 



\n) denotes an exact eigenstate of the numerically solved 
Schrodinger equation for a finite-size disordered system 
with periodic boundary conditions: {Ho + V)'4'n — En'4'm 
where V{v) — uq Y^^' ~ fo'" ^he short-range dis- 
order model we consider. The scattering locations R; 
and potential signs of uq are random. The Schrodinger 
equation has been solved using a large momentum-space 
cutoff fc* « 7 • 10^ cm^^ which corresponds to the en- 
ergy scale at which the split-off bands of bilayer graphene 
become relevant and our two-band model is no longer 
applicable 3^ 

The pseudospin-momentum coupling (the effect in 
which we are mainly interested here) always occurs in 
graphene whichever disorder potential is assumed. The 
model considered here should be seen as a generic one 
where delta-correlated scattering potential is chosen just 
for the sake of simplicity even though the short range 
disorder mixes states in different valleys. The interval- 
ley scattering appears to be irrelevant to any other type 
of disorder (charged impurities, ripples) and is there- 
fore neglected here. Note that the Thomas-Fermi theory 
has been recently employed^^ to calculate the electronic 
structure in the presence of the disorder potential due 
to charge impurities in gapped bilayer graphene. The 
theory is quasiclassical and does not include the quan- 
tum mechanical entanglement considered here. Most 
important is that the amplitude of the screened disor- 
der potential fluctuations must be of the order of the 
gap size A in order to explain the difference between 
the spectral band gap and the experimentally extracted 
transport gap. Here, quite an opposite situation is con- 
sidered: The scatterer strength uq and concentration 
Us = Ns/L? are chosen to be small enough [uqUs < A) 
to preclude the percolative regimei^ and substantial band 
gap renormalization.^- Such careful choice of disorder pa- 
rameters makes it possible to observe the pseudospin co- 
herence effects. 

The Kubo conductivity ([3]) vanishes at 77 — >■ as well 
as at ^ cxd. As one can see in fig. [3l^inset), there is 
an intermediate region near (jimL'^) / (itrh^) = 10 where 
the conductivity is not too sensitive to ry. It is natural 
to work in this region to estimate the conductivity at 
a given system size L. The length L is chosen to be so 
large that the conductivity curves don't change too much 
with further increasing of L. Fig. |H inset) shows the con- 
ductivity curves for different sample sizes starting from 
jL — 0.45 X 10^^ cm with 0.45 x 10^^ cm step. One can 
see that the difference in the conductivity behavior for 
the lengths |L = 1.35 x 10"^ cm and L = 1.8 x 10"^ cm 
becomes rather small, thus, the latter is chosen to be the 
typical sample size which allows the scaling with L. The 
typical scatterer number is a few hundreds for this 
L. The momentum cut-off k* and L fix the Hamiltonian 
matrix dimension at 3362 x 3362. 

The effective mass m for carriers in bilayer graphene 
has been predicted to be equal to 0.054mcr^ with mp be- 
ing the bare electron mass. The latest measurements on 
suspended bilayer graphene^S indicate the effective mass 
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of 0.028mo. In contrast, the charge carriers in bilayer 
graphene on substrate demonstrate larger effective mass 
about O.Oim^ which turns out to be slightly different 
for electrons and holes. Since the model considered here 
assumes the electron-hole symmetry and the effects pre- 
dicted below do not rely on the particular effective mass 
value such precision seems to be excessive for our pur- 
poses and we choose just the round number m = 0.05mo. 

The zero-temperature conductivity curves depicted in 
figs. I3I4I are smoothed by averaging over an energy in- 
terval containing 10-100 levels, over boundary condi- 
tions, and over several disorder potential realizations.— 
The finite-temperature conductivity demonstrates much 
weaker fluctuations, thus, the results shown in figs. I5I6I 
are averaged just over a few disorder realizations. 



IV. RESULTS 

As one can see from fig. |3l the conductivity does not 
vanish even though the chemical potential /x gets below 
the bottom of the conduction band and the tempera- 
ture is zero. Moreover, the subgap conductivity increases 
with disorder (cf. ref.™). This peculiar behavior can 
be understood in terms of the disorder-dependent quasi- 
particle life-time t and pseudospin decoherence time2^ 
Tdc = fi/2Ek with Ek being the characteristic particle 
energy which equals to either /i or A/2 whichever is 
larger. On the one hand, the interband entanglement 
is obviously weaker for larger energies and stronger for 
smaller band gap sizes. On the other hand, the evanes- 
cent components in the interband entangled states be- 
come more important at shorter distances and count in 
favor of strong disorder. As consequence, the subgap 
pseudospin-coherent conductivity contribution increases 
with Tdc/T — the effect we actually observe in figs. I3|4I 
The upper limit for quasiparticle life-time r (which is 
the same as the momentum relaxation time in presence 
of the short-range disorder potential) can be estimated 
using the Fermi golden-rule at ^ ^ J7 as r ~ 3 • 10" s 
corresponding to the mobility lO'^cm^/Vs for curve (b). 

Looking at the plots in fig. |3] one might still think 
that it is the impurity density of states, rather than the 
pseudospin-momentum coupling, that is responsible for 
finite subgap conductivity. In order to clarify this is- 
sue let us compare the pseudospin-momentum coupled 
model ([U with the decoupled one ([2]). The two mod- 
els have the same density of states but the the eigen- 
state spinors do not depend on the direction of particle's 
motion in Here, either conduction or valence band 
eigenstate once created can propagate through the dis- 
ordered sample without changing its pseudospin orienta- 
tion even though the direction of motion is altered after 
each scattering event, as shown in fig. O The interband 
entangled states do not occur here and the conductivity 
vanishes as soon as the chemical potential reaches the 
bottom of the conduction band, see dotted lines in fig. S) 
In contrast, gapped bilayer graphene demonstrates a sub- 
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FIG. 5: Thermally activated conductivity at different band 
gap size A for gapped bilayer graphene (solid lines) and de- 
coupled band intrinsic semiconductor (dotted lines) described 
by eqs. and ((2| respectively. The chemical potential is 
zero, i. e. it is placed exactly in the middle of the gap. 
The subgap conductivity increases slower with the tempera- 
ture within the decoupled band model. The inset shows the 
competition between the temperature-dependent pseudospin 
decoherence and thermal activation of the electron-hole pairs 
resulting in the non-monotonic temperature dependence of 
graphene's conductivity at smaller band gap. Besides the 
band gap size shown in the plot, all other parameters are the 
same as in fig. [3] for curve (b). 



stantial subgap conductivity at the same parameters. 

Thus, to observe the substantial subgap conductivity 
(i) the pseudospin must be coupled with the particle mo- 
mentum to create the interband entangled states in dis- 
ordered samples and (ii) the system must be pseudospin- 
coherent, i. e. r/rdc must be smaller than one. Note 
that Tdc = (for /i > A/2) decreases with increasing 

fi making the two conductivities in fig. 2] indistinguish- 
able at higher carrier concentrations. On the other hand 
the quasiparticle life-time r is longer in higher mobility 
samples with less impurities and/or lighter carriers that 
requires longer Tdc to fulfill the pseudospin coherence cri- 
teria. 

The difference between bilayer graphene described by 
eq. ([H and its rival with decoupled bands ([2]) at best can 
be seen in the thermally activated conductivity. The cal- 
culations can also be considered as a simulation of the 
charge transport in a field effect transistor turned to the 
"off" state when the chemical potential is placed exactly 
in the middle of the band gap hampering both electron 
and hole transport at low temperatures. As one can see in 
fig. [5l the pseudospin-coupled carriers can be excited eas- 
ier than the decoupled ones. The difference between con- 
ductivities in these two cases becomes essential at room 
temperatures. Note that if T ^ A, then the pseudospin- 
incoherent conductivity can be well described by the clas- 
sical formula ^ exp(— ^) indicating that the thermally 
activated conductivity always increases with tempera- 
ture. In contrast, the subgap pseudospin-coherent con- 
ductivity decreases as soon as T becomes comparable 
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FIG. 6: This figure demonstrates the non-monotonic behavior 
of thermally activated conductivity for bilayer graphene at 
the intermediate band gap size A = 43.5 meV for different 
disorder concentrations ria. The chemical potential is zero, 
and disorder parameters for each curve are the same as in 

fig. El 

with A/2 substituting the latter in the expression for Tdc 
and breaking down the pseudospin-coherence. The com- 
petition between these two mechanisms can result in the 
non-monotonic temperature dependence of graphene's 
conductivity, see fig. [5j inset). Note that if T A, then 
both conductivity curves coincide. (This regime is not 
shown in figure.) 

The non-monotonic conductivity behavior is robust 
under moderate change of the disorder strength, see 
fig. ini However, as it was mentioned before, the dis- 
order strength n^uo must always be smaller than the 
band gap size in order to preclude the influence of midgap 
states. The bilayer samples must therefore be relatively 
clean to observe the non-monotonic conductivity behav- 
ior predicted here. The necessary quality can probably 
be achieved in graphene on boron nitride. It is also 
important that the phonons, which are not considered 
here at all, might spoil the effect. The phonon resistivity 
contribution in bilayer graphene is dominated by fiexural 
phonons and rapidly increases with temperature. « The 



fiexural phonons can be again suppressed in graphene 
sandwiched between boron nitride layers i^"— 



V. CONCLUSION 

To conclude, there is a fundamental obstacle which 
limits the functionality of the field effect transistor based 
on gapped bilayer graphene. The physical mechanism re- 
sponsible for that is intimately linked to the pseudospin- 
momentum coupling which leads to the instantaneous 
generation of the interband entangled states in the pres- 
ence of disorder. It makes higher "leakage" current in 
the "off" state and therefore limits the possible on/off 
ratio by lower values as compared to those in conven- 
tional semiconductor devices with the same mobility and 
band gap size. In contrast to the "leakage" mechanisms 
considered before r^ii^ii^ the interband entanglement de- 
scribed here is unavoidable unless the very crystal lattice 
is broken. Moreover, in contrast to the universal subgap 
conductivity observed in the topological insulators the 
subgap conductivity in bilayer graphene turns out to be 
sensitive to the band gap size and disorder strength. The 
non-monotonic conductivity vs. temperature dependence 
predicted here can be seen as a signature of the pseu- 
dospin precession responsible for the difference between 
the transport and spectral gaps. The effect can probably 
be observed in doubly gated bilayer graphene sandwiched 
between boron nitride layers where the charge inhomo- 
geneity and fiexural phonon conductivity contributions 
are substantially reduced. 
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